Mon. Not. R. Astron. Soc. OOP. [TlfTOl (20111 Printed 16 September 2011 (MN WI^ style file v2.2) 



Extension of an Exponential Light Curve GRB Pulse 
Model Across Energy Bands 



o 
in 



o 



> 

m 
m 

O 



Robert J. Nemiroff^* 

^Department of Physics, Michigan Technological University, Houghton, MI 49931 



ABSTRACT 

A simple mathematical model of GRB pulses in time, suggested in lNorris et alj (|2005[ ). 
is extended across energy. For a class of isolated pulses, two of those parameters 
appear effectively independent of energy. Specifically, statistical fits indicate that pulse 
amplitude A and pulse width r are energy dependent, while pulse start time and pulse 
shape are effectively energ y independent. T hese results bolster the Pulse Start and 
Pulse Scale conjectures of iNemirofj (j2000l ) and add a new Pulse Shape conjecture 
which states that a class of pulses all have the same shape. The simple resulting pulse 
counts model is P{t,E) = A{E) exp{—t/T{E) — T{E)/t), where t is the time since 
the start of the pulse. This pulse model is found to be a n acceptable statistical fit 
to many of the fluent separable BATSE pulses listed in iNorris et al.l (|2005[ ). Even 
without theoretical interpretation, this cross-energy extension may be immediately 
useful for fltting prompt emission from GRB pulses across energy channels with a 
minimal number of free parameters. 
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1 INTRODUCTION 

Gamma Ray Burst (GRB) prompt emission appears typi- 
cally to be co mposed of d i stinct emission epis odes known 
as "pulses" (|Desail Il98ll : iNorris et al] 1 19961 ). Although 
pulses typically overlap in time, a fraction of GRBs fea- 
ture a pulse bright enough, long enough, and separate 
enough to be analyzed by itself. Previously, several authors 
have sugges ted relatively simp le analytic forms for GRB 
pulses (e.g. iNorris et ail (|l996l ): iRvde fc SvenssonI (|2002h : 
INorris et al.l ( 20051 ') 'I. and several conjectures have been 



given and tested involving GRB pulse coherence (|Nemirofll 
I2OOOI : iHakkila fc Nemir"oa (20091 ). 

The INorris et al.l (| 19961 ) and INorris et al.1 (|2005l ) math- 
ematical forms fit GRB pulse light curv es to an exponen- 
tial ri se and exponential decay, while the lRvde fc SvenssonI 
l|2002h form fits GRB pu lse light c urves to a power law 
rise and power law decay. iNemiro^ (|201ll ) fit a Planckian 
functional form that has an expone ntial rise and a power- 
law decay ISchaefer fc DvsonI (|l996h showed that 10 GRB 
pulses detected early in the Burst and Transient Source Ex- 
periment's (BATSE's) mission, on board the now-defunct 
Compton Gamma Ray Observatory, were only marginally 
well fit to an exponential decay, while a power-law decay 
sometimes fits better. 

Her e a simple analytica l pulse models suggested previ- 
ously by INorris et al.l (|2005l ) is analyzed in greater detail. In 
particular, the pulse model is reformulated mathematically 
into three component parameters; one that scales linearly 



in time, one that scales linearly in intensity, and one that 
solely determines the pulse model shape. These parameters 
easily scale with energies at which the pulse is observed. 

All previously published pulse fitting schemes at- 
tempted to fit a pulse only at a single energy or in a single 
energy band. Fits to the same pulse at another energy are 
typically started fresh, with all the free parameters again be- 
ing determined from scratch. To date, no published system 
uses information from a pulse at one energy to fit the same 
pulse at another energy. Deconvolving complicated GRBs 
into pulses , however, can be a computationally expensive 
procedure (|Hakkila et al.|[2008l ). The problem is not just an 
inefficient use of computer time - it affects fitting accuracy 
as well - information gained from fitting the pulse in a bright 
energy channel is not being used to formally constrain the 
fit of a dim en ergy channel, altho ugh it may be used as a 
starting point (jHakkila et al.ll2008l ). 



In Section 2 the mathematical [Norris et aL 1 (|2005h pulse 
model is reviewed and reformulated in terms of more eas- 
ily scalable param eters. Section 3 shows how the scalable 
INorris et al.l (|2005l ) pulse is extendable across energy chan- 
nels. Section 4 shows fits for the energy generalized pulse 
mod el to several bright individual BATSE pulses selected 
from INorris et al] (|2005l ). Section 5 details a search for cor- 
relations between scalable parameters found in the fits, while 
in Section 6 some discussions and conclusions are given. 
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2 THE NORRIS PULSE MODEL 

2.1 The Norris Paramatrization 

In lNorris et"aD (j2005h . the following pulse model was fit to 
numerous GRB pulses and discussed: 

P{t) = Ae-^'-'^-^'''\ (1) 

where P is the count rate of the pulse in counts per second, 
t is time during the GRB, A is the amplitude of the pulse 
also in counts per second, Tt is a temporal factor that scales 
the rise of the pulse, and Td is a temporal factor that scales 
the decline of the pulse. Here the pulse start time to has 
implicitly been set to zero. Although the first author of this 
paper, J. P. Norris, has been involved in the creation and 
testing of other pulse models, here this puls e model will 
be ref erred to as the "Norris pulse model". In lNorris et al.l 
(|2005l ). pulses were fit separately to the Norris pulse model 
in every energy band. 

As with any peaked function, the time of peak counts 
is found simply by finding when dP/dt = 0. For the Norris 
pulse model this occurs when t = tp^ak ~ v'tvt^. At this 
time, the maximum count rate of the pulse is 

Ppeak = P{tpeak} = A e-"^/^. (2) 

The pulse shape can be uniquely determined by also find- 
ing P at another time, here chosen to be at twice the time 
between the peak and the start of the pulse, so that 

P{2tp,ak) = A e-(5/2)^/W^. (3) 

2.2 A Different Parametrization of the Norris 
Pulse Model 

The same pulse shape can be rewritten in terms of other vari- 
ables with the goal of giving greater insight into the cross- 
energy scaling potentially inherent in GRB pulses. Specifi- 
cally, Eq.fT]) can be rewritten as 

P(f) = Ae-«'*/"+"/". (4) 

Since r — ^/ryrd and ^ = a/ Tr/ra, it is clear that this pulse 
description involves just a direct variable substitution and so 
is mathematically identical to the previously defined Norris 
Pulse Model. This pulse light curve parametrization appears 
more easily scalable than the original parametrization, a fea- 
ture that will be discussed in some detail below. 

In this parametrization, the peak of the pulse occurs at 
tpeak = T where 

Ppeak ~ P{ipeak) ~ Ac ^, (5) 

so that 

5 = -{l/2)ln{P{tp,ak)/A). (6) 

Note that this does not uniquely specify ^ since the value of 
A remains unknown, and the value of tpeak is known only 
relative to the start time of the pulse: to- 
At twice the peak time, however. 

Pa = P(2fpe,fc) = Ae-^«/^ (7) 

so that Ppeak /P2 = c^^^'^ . Siuce this ratio eliminates A, it 
is now possible to show that 

C = -2 HPpc^ak/P2), (8) 
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Figure 1. A plot of count rate versus time for theoretical Norris 
pulse model for different amplitudes A. Here the pulses sliare 
all other defining parameters, specifically § = 1, r = 100, and 
to = 100 sec. An artificial background level of 100 counts per 
second was input. Note how changing A only scales the count 
rate axis (here the y-axis), while leaving all other aspects of the 
pulse unchanged. 



which does uniquely specify ^ in terms of the relative value 
intensity of the pulse peak to the intensity at twice the time 
of the peak. 

This scalable parametrization allows the deconvolution 
of the pulse shape into three parameters that have separa- 
ble and effectively orthogonal meanings. First, scaling A will 
stretch the counts ("y-axis") of the pulse only, leaving the 
timing ("x-axis") and the inherent shape of the pulse un- 
touched. This is shown graphically in Figure ([T}. Next, scal- 
ing T will stretch only the time axis of the pulse, leaving the 
amplitude and the inherent shape of the pulse untouched. 
This is shown graphically in Figure ([2]). Note that this time 
stretching occurs while holding the pulse start time to fixed. 

Last, in this deconvolution, the ^ parameter alone de- 
fines the scale-independent shape of the pulse. Only when ^ 
is changed, for example, can the asymmetry of the pulse be 
changed. This is shown graphically in Figure ([Sjl. Changing 
the shape parameter ^, however, will also change the ampli- 
tude of the pulse. The only remaining free parameter is to, 
the start time of the pulse. 

Note that in this parametrization, r, the time of the 
peak relative to the start time of the pulse, can b e used as a 
proxy for width para meter w used InlNorris et al.l l|2005h . As 
shown in Figure (4) of iNorris et al] (|2005l l. there is excellent 
agreement between the two. 

Since tpeak = r, iL straightforward method for estimat- 
ing T is between the measured time values of the pulse 
peak and to. Although measuring the peak time is rela- 
tively straightforward, measure the start time is more com- 
plex. However, since most pulses rise relat i vely r apidly com- 
pared to their decay, iHakkila fc Nemirofj (|2009t ) , for exam- 
ple, found it possible to estimate iovalues to an accuracy 
below 0.4 seconds from a statistical fit of the Norris pulse 
form. 
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Figure 2. A plot of count rate versus time for theoretical Norris 
pulse model for different time scaling factors r. Here the pulses 
share all other defining parameters, specifically ^ = 1, A = 60, 
and to = 100 sec. An artificial background level of 100 counts 
per second was input. Note how changing t only scales the time 
axis (here the x-axis), while leaving all other aspects of the pulse 
unchanged. 
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Figure 3. A plot of count rate versus time for theoretical Nor- 
ris pulse model for different pulse shape parameters ^. Here the 
pulses share all other defining parameters, specifically A = 60, 
T = 100, and to = 100 sec. An artificial background level of 100 
counts per second was input. Note how changing ^ fundamentally 
changes the shape of the pulse. Note that lower values of ^ cre- 
ate more asymmetric pulses, while higher values of ^ create more 
symmetric pulses. 



3 EXTENSION OF NORRIS PULSE MODEL 
TO OTHER ENERGIES 

The scalable parametrization of the Norris pulse model can 
be easily generalized to other energies, so long as the pulse 
light curve is fit acceptably at those energies. Generally, Eq. 
(H} can be expanded so that 



(E) ' (t-to(E))y 



(9) 



where it is now assumed that A, r, 5, and to are all functions 
of energy E. 

This generalization would be particularly useful if any 



of the functions A{E), 5(e), t{E), or to{E) were known since 
that would reduce the number of free parameters needed to 
fit a pulse. Such a function might come from theory one 
day, but presently one can postulate a few empirical rules 
bounding these functions from the inspection of numerous 
GRB light curves. One such empirical hypothesis is that 
the values of these functions are smoothly varying over E, 
so that the knowledge of a functional value at one energy is 
likely similar to a value at a close energy. Past this, however, 
only two concrete relations have been hypothesized to relate 
GRB pulse light curves be tween energies: the Pulse Start 
Conjecture (lNemiroflll2000l ) and the Pulse Scale Conjecture 
(|Nemirofj|2000l V 

The Pulse Start Conjecture posits that every pulse has 
the same start time at every energy, the mathematical state- 
ment of which is that to{E) = to, a, constant that is not a 
function of en ergy. The Pulse Start Conjecture was posed 
and tested in iNemiro^ (|2000l ) for the most fiuent pulse 
in BATSE trigger 2193 (GRB 930214c) for 16 MER en- 
ergy channels and found to be an acceptable statistical fit 
for all of them. Further, the Pulse Start Conjecture was 
tested in Ncmirofi (200Qi) using BATSE four-energy chan- 
nel data for the most fluent pulses in four other BATSE 
long GRBs (GRB 930612a, GRB 940529b, GRB 941031b, 
and GRB 970825, corresponding to BATSE trigger num- 
bers 2387, 3003, 3267, and 6346, respectively) and found to 
be statisticall y consistent. 

However. iNorris et al.l (|2005l 'l measured to values for the 
brightest pulse in four energy channels for BATSE GRB 
2193 while fitting the pulse. The most discrepant fit was 
found to be start times that differed by 2.0 a from each 
other. A statistically significant difference was recorded for 
BATSE trigger 2387 between energy channels 1 and 2, and 
channels 1 and 3, where the iNorris et al.1 (|2005l ) to values 
differed by 4.7 a and 8.2 a respectively. However, these dif- 
ferences could be due to a single displacement i n the channel 
one st art time for this pulse. More generally, INorris et al.l 
(|2005h measured start values for numerous bright pulses in 
BATSE GRBs and recorded data which include cases that 
appear to follow the Pulse Start Conjecture, and some cases 
that don't, assuming that the error in to measured is gaus- 
sian. 

Most recently, iHakkila fc Nemiroflj (|2009l ') tested 199 
pulses found in 75 long GRBs and and found that the Pulse 
Start Conjecture to be generally plausible, with a few errant 
cases possibly attributable to secondary pulses or low signal 
data. 

Given that the Pulse Start Conjecture holds, then Eq. 
((9)1 can be simplified to 



«(^5)(7I1 



(10) 



where t is now considered measured relative to an energy 
independent to- 

A separate conjecture, the Pulse Scale Conjecture, 
posits that GRB pulses differ between energies only by scale 
factors in time and flux (Ncmiroff .2000i ) . In this paradigm, 
it is 5 that is constant and hence no longer a function of 
energy. Note that the Pulse Scale Conjecture allows differ- 
ent pulses to have different shapes. Given that Pulse Scale 
Conjecture holds, then Eq. (|9]) can be further simplified to 

P(t,i5) = .l(i?)e-«"/^(^)+^(^)''*\ (11) 
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The remaining energy dependant terms are now A(E) and 
r{E), with A{E) now usefully being the "y-axis" flux scale 
factor and r{E) being the othogonal"x-axis" temporal scale 
factor. 

A yet further conjecture ventured here, the Pulse Shape 
Conjecture, posits that a group of GRB pulses all have the 
same shape. In this paradigm, ^ is not only constant in a 
pulse across energies, an entire group of pulses all are hy- 
pothesized to have the same shape: ^ = 1. Given that the 
Pulse Shape Conjecture holds, then Eq. (|9]) can be further 
simplified to 

P(t,i5) = A(i5)e-(*/^(^)+^(^)/'\ (12) 

It is suspected that this functional form is not unique in that 
other mathematical functions can also adequately fit all of 
the pulse conjectures and describe the same group of GRB 
pulses. 

In the above cross-energy pulse descriptions, the mean- 
ing of A{E) is interesting as it appears to be a type of spec- 
trum for a pulse that involves the fiux from the entire pulse. 
Note that A{E) will be different than a spectrum taken at a 
specific time for of the pulse, say at the time of peak energy 
of one of the energy channels. Interestingly, A{E) can be fit 
from incomplete pulse information, for example computable 
when only a fraction of a single pulse is recoverable across 
energy bands. 

Previously, A{E) has only been determined for a single 
pulse in a sing l e GR B: the main pulse in BATSE trigger 
2193 (|Nemirofj[2000l ). There, A{E) was fit to 16 BATSE 
MER channels and was seen to peak at an energy of about 
150 keV. Note that A{E) is intrinsically a type of "photon 
count spectrum" , as P{t) is a measurement of photon counts 
rather than energy fiux. 

The meaning of t{E) appears to be the pulse dura- 
tion as a function of energy. Note that t{E) has (also) only 
been determined for single pulse in a single GRB, i n fact the 
same pulse (BATSE trigger 2193) as with A{E) (|Nemirofj 
|2000D . There, t{E) appears approximated as a power law 
where r{E) oc _E~°'^. Informal observation of numerous 
GRB pulses indicates that t{E) usually decreases mono- 
tonically with increasing energy. Note that as with A{E), 
r{E) may also be fit from incomplete pulse information, for 
example computable when only a faction of a single pulse is 
recoverable across energy bands. 

Note that the total fiuence of a Norris-model pulse is 
directly proportional to A{E) times t{E). Since A{E) has 
a clear peak, pulses might be seen as envelopes effectively 
spanning only certain times and energies. 



4 FITTING GRB PULSES WITH 

CONSTRAINTS ACROSS ENERGY 
CHANNELS 

In this section, an attempt is made to determine whether 
the Norris pulse model can be constrained by the Pulse 
Start, Pulse Scale, and Pulse Shape conjectures and still ad- 
equately fit isolated GRB pulses adequately across BATSE 
energy bands. Note that invoking these conjectures provides 
many additional constraints, so that many fewer free param- 
eters will be available to fit the pulses across the BATSE 
energy bands. As a consequence, one might expect fits to 



be less good than fits without these constraints. Therefore, 
what is being tested is not whether such constraints pro- 
vide better fits - but whether these more tightly constrained 
models can create acceptable fits at all. 

BATSE data was used for several reasons. First, BATSE 
incorporated the largest set of GRB detectors ever ffown, 
and so recorded relatively high count rates for GRBs. Next, 
BATSE, so far, was one of the longest running GRB detec- 
tors ever ffown, and so over its decade of operation had time 
for wait for relatively bright GRBs to occur. Next, BATSE 
data has now been around for over two decades and so is rel- 
atively well understood. Next, BATSE data is in the public 
domain and easily availab le over the web in s imple formats, 
for example ASCII. Last, iNorris et al.l (|2005l ) created their 
bright isolated pulse catalog from BATSE data, which is a 
catalog with useful data from which this study can draw and 
compare. 

In general, four channel BATSE LAD data were used. 
The four channels used were channel 1 (20 - 50 KeV), chan- 
nel 2 (50 - 100 KeV), channel 3 (100 - 300 KeV), and channel 
4 (300 - 1000 MeV). In many GRBs, the counts in channel 
4 are so low that the burst may not be easily discernible,. 
Therefore, only the fits for the first three BATSE LAD en- 
ergy channels are reported in Table 1. The BATSE LAD 
data type with 64-ms time resolution was used, as it had 
the highest temporal resolution that covered greatest tem- 
poral range in and around pulses. 

The procedure for determining which GRB pulses were 
selected for fitting to the scalabale Norris pulse f u nction was 
as follows. All of the GRBs in the iNorris et all (|2005h cat- 
alog were c onsidered initially. These pulses were considered 
initially bv INorris et al.l (|2005l ) because they were isolated, 
of relatively long duration, were relatively ffuent, and had 
relatively long lag times between cross correlations between 
BATSE energy channels. A visual inspection was then done 
excluding some GRBs that were subjectively deemed too 
convolved with other pulses to yield a meaningful result. 
Admittedly, this procedure might bias the results so that 
only pulses that might fit are selected, but this was con- 
sidered acceptable since a key premise being tested is that 
some GRB pulses can be described by the energy-extended 
Norris pulse form, not all of them. 

Once a GRB and its dominant pulse were selected, a 
range of times inside the GRB when the pulse being consid- 
ered was clearly uncontaminated by surrounding pulses were 
recorded. The pulse light curve was then rebinned in time 
so that the pulse rise time occurred between about 10 to 20 
time bins. This allowed the pulse to be temporally resolved 
while retaining a relatively large signal to noise per time 
bin. Two background intervals, typically lasting 100 seconds, 
were selected well before the pulse and after the pulse and 
fit with a single third degree polynomial. A possible range 
of reasonable start times, to, were estimated based on visual 
inspection of the pulse light curves in the four BATSE LAD 
energy channels. 

To test for goodness of fit, candidate to values in this 
range were assumed, in turn. For each to assumed, the best 
fit T and A values were recorded that minimized the per 
degree of freedom {v, the number of time bins) for each 
energy channel. The value of to that allowed the best fit A 
and r values were recorded, as well as error estimates for 
these values. 
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Table 1: Temporal Fit Parameters 
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990220 


7403 


100 - 


300 


-3.6 


375 


6.2 


1.67 
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Table 1 Continued 



GRB 


BATSE Trigger 


Energy Range 


to 


A 


T 




990528 


7588 


20 - 50 


-1.3 ± 0.4 


267 ± 19 


4.93 ± 0.23 


0.99 


990528 


7588 


50 - 100 


-1.3 ± 0.4 


350 ± 25 


4.11 ± 0.19 


1.26 


990528 


7588 


100 - 300 


-1.3 ± 0.4 


267 ± 39 


3.13 ± 0.29 


1.13 


990712B 


7648 


20 - 50 


-2.2 ± 0.6 


144 ± 21 


7.78 ± 1.56 


1.05 


990712B 


7648 


50 - 100 


-2.2 ± 0.6 


189 ± 28 


7.10 ± 1.05 


0.90 


990712B 


7648 


100 - 300 


-2.2 ± 0.6 


249 ± 18 


5.40 ± 0.24 


1.10 


000323 


8049 


20 - 50 


-3.4 


249 


14.73 


1.17 


000323 


8049 


50 - 100 


-3.4 


402 


11.21 


1.44 


000323 


8049 


100 - 300 


-3.4 


430 


8.15 


2.00 



There are many factors that can affect the statistics of 
the fit so as to render the computed goodness of fit values 
more indicative than numerically definitive. These include 
the fineness of the time binning, the precision of the back- 
ground fit, the temporal range of pulse bins fit, and the affect 
of secondary independent pulses, fn general, breaking the 
pulse up into more time bins tended to lower ^ jv^ although 
the increase in v then tended to offset significant changes in 
actual goodness of fit. Although pulse fits always started at 
to, extending the temporal fit range of the pulse decay fur- 
ther into the background tended to drive jv down since 
the more slowly changing background was typically easier 
to fit. The pulse was deemed to have ended, in terms of the 
statistical fitting procedure, when the pulse was either a few 
percent above background or a second contaminating pulse 
became clearly apparent. Given these boundaries and limi- 
tations, a pulse fit was considered descriptively useful when 
X^/i' was computed to be 1.5 or lower over the fitted pulse 
range. 

The results of the fitting process are shown in Figure 
(|4]) with the best fitting parameters listed in Table 1. For 
a GRB's considered pulse, a best fit was declared across 
energy channels when the sum of the squares of 'x^ jv over 
all of the energy channels was at a minimum, fn Figure 
(|4)|, fight curves for the GRB are given for the lowest three 
BATSE energy bands, with the lowest energy given at the 
top. The raw BATSE data is shown in histogram form by the 
fluctuating solid line. Time is depicted relative to the start 
time of the dominant pulse shown, with a vertical dashed 
line making this start time to more clearly. Superposed on 
the data are the best fit plots for the Norris pulse model 
expanded across energy bands as described above. 

In Table 1, the first column lists the GRB date desig- 
nation, while the BATSE trigger number is given in column 
2, and the BATSE energy range being fit is listed column 3. 
The best fit start time to relative to the trigger time for the 
dominant pulse in the GRB is given in column 4, along with 
an error estimate. For this start time, the best fit parameters 
for A(E) and t{E) are given in columns 5 and 6 respectively, 
along with their error estimates. The jv goodness of fits 
values are listed in column 7. 

Error estimates are only given for GRBs where thex^/i^ 
was below 1.5 for every energy channeL For GRB pulses 
where the fit was deemed "not good," meaning here that at 
least one energy channel could not be fit with the minimal 
Norris pulse function to a ^ jv better than 1.5, then only 



the best fit to, A, and r parameters are listed. Where a 
good fit was found, the to. A, and t values, as well as their 
errors, are listed. The errors indicate the limiting value of 
the parameter that could create a good fit, given the best 
fit listed for the other parameters. 

For BATSE 6598, a better fit was obtained for ^ > 1. 
Specifically ^ = 2 provided a better fit in almost all en- 
ergy channels. BATSE trigger 6625 appears upon inspec- 
tion to contain an early pulse that is biasing the statistics. 
In BATSE trigger 8049, only the first pulse was fit. This 
is an example of fitting a pulse based on incomplete infor- 
mation, which this technique allows. Note that for BATSE 
8049, the second pulse was not used since one would first 
have to model the first pulse and subtract it, which was con- 
sidered beyond the scope of this work. Given the Pulse Start 
Conjecture, though, the second pulse is hypothesized to be 
completely contained after a time that can be estimated by 
inspection of the light curve. 



5 A SEARCH FOR CORRELATIONS 
BETWEEN SCALABLE PULSE 
PARAMETERS 

From inspection of Table 1 and the plots in Figure ), it is 
clear that a group of GRBs exist that contain a pulse that 
statistically satisfies all three of the above stated pulse scale 
conjectures simultaneously. From the GRBs tested here, this 
group includes BATSE trigger numbers 2197, 2665, 3256, 
5387, 5415, 6147, 7087, 7588, and 7648. The dominant pulse 
of these well-fit GRBs will be called an Energy Scalable 
Pulse (ESP). 

Inspection of the A values in Table 1 shows that the 
full pulse spectra for ESPs can be quite diverse. Typically, a 
spectrum will peak in one of middle BATSE energy bands, 
hence showing an Ep^ak value common to the range of Ep^ak 
values in published GRB spectra. The existence of Ep^ak in 
this range might not be coincidence and might be related 
to the trigger criteria of BATSE GRBs. The low energy res- 
olution of the spectra make it difficult to determine if the 
full pulse spectra follow a Band functional form or any other 
spectral form common to other types of GRBs. 

A plot depicting how the temporal stretch factor r 
changes with energy for the ESP sample is shown in Figure 
©. Each GRB is labelled by its BATSE trigger number. 
Note that r is usually consistent with a power law in energy 
to within the uncertainty. There is at least one exception to 
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Figure 4. Plots of GRB light curves. 
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Figure 4. Continued 
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Figure 4. Continued 



o 



31.6 



10.0 



3.2 



1.0 



10 




32 



100 
E (KeV) 



316 



1000 



Figure 5. A plot of pulse stretehing faetor r versus pulse energy 
E for selected pulses where where the Norris pulse model approxi- 
mations provided reasonably good data fits to BATSE data. Error 
bars are included, where the error in E was just the width of the 
BATSE energy channel. Each point is labeled by its BATSE trig- 
ger number. Note that for most GRB pulses fit, t may be related 
to E by a power law. 



atively large difference between the duration of the pulse 
at high and low eneries. In contrast, the dominant pulse in 
BATSE trigger 7087 appears weU fit by r cx E-°-^'^ , yield- 
ing a relatively small difference between the duration of the 
pulse at high and low energies. 

Note that the difference between the r factors of two 
energy channels is conceptually similar to the lag between 
these two energy channels. It is therefore possible that the 
power law exponent between r and E might act similar to 
pulse lag and be a measure of intrinsic brightness of the 
pulse and hence the GRB. Since r is the time between the 
pulse peak and the pulse start, the difference between two 
rs of the same pulse at different energies is just the dif- 
ference between the peak times of the pulse at those ener- 
gies. Although the lag is formally computed by noting the 
time of the maximum in the cross correlation between the 
light curves of two energy channels, this cross correlation 
would have it's greatest instantaneous power when the two 
peaks are aligned. Therefore lag is strongly pulled toward 
the difference in peak times. Therefore, one would expect 
lag23 ~ (ts - n). 



6 DISCUSSION AND CONCLUSIONS 



this rule, in that BATSE trigger 2665 does not appear to 
have a power law relation between r and E. 

Furthermore, the exponent of the power law relation- 
ship between E and r does not appear to be consistent be- 
tween GRBs. For example, the dominant pulse in BATSE 
trigger 2197 appears well fit by r oc yielding a rel- 



it is evident from inspection that GRB pulses have simi- 
larities at different energies. Here clear rules are postulated 
for transferring specific pieces of information between en- 
ergies. Specific analyses presented above indicate that at 
least one class of GRB pulses, here dubbed "energy scalable 
pulses" (ESPs), can be described by a simple, cross-energy 
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mathematical form that has a well-defined start time and 
a common light curve shape that difi'ers between energies 
by only scale factors is time and brightness. The particu- 
lar li£lit_cuTwslM£e ex here is the simple one given 
in iNorris et all (|2005l ): P{t) = yle-(*/^+^'''' where t is time 
since an energy independent pulse start to, and (only) A and 
r are functions of ener gy. Other simpl e light curves descrip- 
tions are also possible ( Nemiroflll201lh . 

Implications of the above results can be broken into 
two categories: computational and physical. Computation- 
ally, ESPs can be adequately fit to real data by exploring a 
reduced parameter space when compared to non-ESPs. For 
example, fitting a non-ESP, or ignoring the ESP nature of 
a pulse, may cause it to be fit by four free parameters per 
pulse per energy channel: to, A, r, and ^. Invoking the ESP 
model halves the number of fit parameters to explore, since 
start time to and pulse shape ^ are postulated to be the 
same at each energy. The time saved computationally may 
exceed a factor of two since the number of nested loops is 
halved once to and ^ have been determined. Furthermore, ^ 
may not even have to be determined when postulated to be 
same for different pulses. 

The physical implications of the existence of ESP pulses 
might be numerous, but only a few aspects will be touched 
on here. First, the scalable nature of complete ESPs be- 
tween energies demands that the rise and fall of the pulse is 
dependent on the same temporal scale factor r, p erhaps in- 
dicat i ng a common energy lib eration mechanism iNemirofj 
I2OOOI: iHakkila fc Preecjl201ll ). This is counter to analyses 



Norris et al.l (|2005h where the pulse rise and fall are con- 



sidered separate, implying that the energy liberation mecha- 
nisms might also be separate. It is still possible, for example, 
that a correlation exists between duration r and pulse shape 
^ for non- ESP pulses jPeng et all I2OI0I : iHakkila fc Preecd 

Next, the common start time to found here for ESPs 
across energy bands further bolsters the case that GRB 
pulses are generated by some physical triggering mechanism 
related to that start time. Whether to marks the time of 
an explosion, a collision, a phase change, or something com- 
pletely different, is not speculated on here. 

It is interesting that the error in this start time is typi- 
cally much less than the duration of the pulse, typically by 
at least one order of magnitude. In the case of the domi- 
nant pulse in BATSE GRB 5415, the ratio of Afo/r was 
recorded to be about 0.055 for energy channel 1, while Ato 
was estimated to be about 0.4 sec. Seemingly trivial re- 
sults like this might be of particular interest because they 
might provide important limits on any dispersive property 
of the universe between us and cosmo logically distant GRBs 
(jNemiroff. Holmes, fc Connelly! [20l3 '). 

Next, the common shape parameter ^ = 1 shared by 
many ESPs may indicate that a single set of physical pa- 
rameters may be operating on these pulses. Physical mech- 
anisms that are capable of generating a shape statistically 
similar to the Norris pulse shape should be explored. 

Why aren't all GRB pulses ESP pulses? It is possible 
that they are, but that the temporal overlap between multi- 
ple pulses is hiding this a t tribut e for most pulses. In many of 
the cases ofH orris et all (|2005l l selected pulses, a dominant 
pulse might be contaminated by at least one much smaller 
amplitude pulse. One possible example of this is BATSE 



trigger 6625, shown in Figure 4. There, most of the late time 
behavior of the pulse appears well fit by the ESP paradigm, 
while a faint pulse may contaminate the start of the rise 
period of the dominant pulse. 

Possibly, whatever processes that create GRB pulses ac- 
tually generate a whole luminosity function of pulses with 
lower amplitudes and lower energies being more common. 
Then, most detected GRBs are such dense conglomerations 
of these faint and weak pulses that it is not possible to sep- 
arate out even one pulse for detailed analysis. The selection 
of the lNorris <^~^ (|2005l ) sample was an attempt to isolate 
single separable pulses, but it is understandable that other 
pulses might have been inadvertently included, in particular 
faint pulses. 

It does appear, however, that some pulses are just not 
good fits to the ESP paradigm. It is interesting to wonder 
if all the conjectures that go into the Norris pulse fit are 
violated, if some of them are violated, or perhaps just one 
goes awry. Some conjectures may also be more true, or better 
approximations, than others. 

From the sample fit, it is speculated here that the pulse 
start conjecture is the strongest approximation, the pulse 
scale conjecture is the next strongest, and the pulse shape 
conjecture is the weakest, in particular between different 
pulses. In other words, it appears that some GRB pulses are 
better fit with ^ 7^ 1, but that this shape may be consistent 
across energy bands for that pulse. Additionally, it appears 
that other pulses just cannot be fit with any ^, being, for 
example, too sharply peaked in the center. 

For well fit ESPs, it is hoped that this work has better 
quantified not only how A and r change with energy, but 
how better known parameters, such as peak flux (propor- 
tional to A) and fiuence (proportional to At) change with 
energy. When a GRB ESP sample with actual redshifts be- 
comes known, it is hoped that this energy-dependent cor- 
relations might one day lead to lower dispersion standard 
measures bolstering the use of GRBs as cosmological probes. 

As with much research, this work might just be the tip 
of a much larger iceberg. Many questions arise that are not 
yet answered, and might take considerable effort and greater 
data to answer. For example, over what energy range does 
the ESP paradigm work? Perhaps GRB data over a wider 
energy band would indicate clear energy cutoffs at both high 
and low energies where pulse shapes do some something 
completely different. 

The author wishes to acknowledge useful discussions 
with Jay P. Norris, Jerry T. Bonnell, Demos Kazanas, 
John Hakkila, Amir Shahmoradi, Daniel Miller, and Justin 
Holmes. 
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